Microcanonical mean-field thermodynamics of self-gravitating and rotating systems 
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We derive the global phase diagram of a self-gravitating iV-body system enclosed in a finite three- 
dimensional spherical volume V as a function of total energy and angular momentum, employing a 
microcanonical mean-field approach. At low angular momenta (i.e. for slowly rotating systems) the 
known collapse from a gas cloud to a single dense cluster is recovered. At high angular momenta, 
instead, rotational symmetry can be spontaneously broken and rotationally asymmetric structures 
(double clusters) appear. 
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The statistical equilibrium properties of systems of 
particles interacting via long-range forces (the so-called 
non-extensive systems) are currently the subject of in- 
tense research, both for their highly non-trivial thermo- 
dynamics (displaying such features as negative heat ca- 
pacities 111 |2|) and for the considerable conceptual and 
technical difficulties they present. It is known that the 
long-range nature of the potential makes the canonical 
ensemble inadequate for describing their statics || |[ |), 
because the usual thermodynamic limit, where (number 
of particles) N — > oo and (volume) V — > oo while in- 
tensive variables are kept fixed, does not exist. A cen- 
tral issue is hence whether phase transitions and other 
conventional statistical phenomena are possible in non- 
extensive systems ||. 

Among non-extensive systems, self-gravitating gases, 
i.e. systems of classical particles subject to mutual grav- 
itation, have deserved the most attention. Their usual 
static description is based on the microcanonical ensem- 
ble Q. In this framework, the key problem is finding 
the most probable equilibrium configuration of a self- 
gravitating gas enclosed in a finite 3-dimensional box of 
volume V as a function of the conserved quantities (inte- 
grals of motion), the simplest (but possibly not the only 
relevant ones) being the total energy E and the total 
angular momentum L. 

Dynamical methods H based on fluid-mechanics tech- 
niques suggest (see e.g. ]S ]) that upon increasing the ratio 
between rotational and gravitational energy, the station- 
ary distribution can change from a single dense cluster to 
a double cluster, and that other structures such as disks 
and rings might appear. 

On the other hand, so far static theories could not re- 
cover the richness of the dynamical picture. Taking the 
total energy as the only control parameter (see |g for 
a review and || [h], [ll| for more recent work and refer- 
ences) after removing the rotational symmetry artificially 
e.g. by constraining the system into a non-spherical box, 
a "collapse" transition has been found where, as the 
energy (temperature) is lowered, the equilibrium configu- 
ration changes from a homogeneous cloud to a dense clus- 
ter lying in an almost void background, with an interme- 
diate "transition" regime characterized by negative spe- 



cific heat. Despite some attempts |Q, a detailed static 
theory embodying angular momentum is lacking. 

In this work, building essentially on ||, we aim at 
bridging the gap between the static and the dynamical 
approaches by analyzing the statics of a self-gravitating 
and rotating gas in a microcanonical setting including 
angular momentum. Within a mean-field approximation, 
we derive an algebraic integral equation for the density 
profiles maximizing the microcanonical entropy and solve 
it numerically as a function of E and L. Along with the 
usual collapse, occurring at low L, we find that for suf- 
ficiently high L the rotational symmetry of the Hamil- 
tonian can be spontaneously broken, giving rise to more 
complex equilibrium distributions, including double clus- 
ters, rings and disks. The global phase diagram of the 
system is presented. We shall concentrate here on the 
equilibrium density profiles, deferring a detailed discus- 
sion of the related (highly non-trivial) thermodynamic 
picture to a more extensive report. 

We consider the A-body system with Hamiltonian 

H N = H N ({r t }, { Pi }) = — J2p1 + $({r *}) (1) 



with $({r,-}) = -Gm 2 Y, i<J k, - fj\ n, p t and m 
denote, respectively, the position, momentum and mass 
of the i-th particle. The system is assumed to be enclosed 
in a spherical volume V (to preserve rotational symme- 
try and ensure angular momentum conservation). The 
crucial quantity to be evaluated is the microcanonical 
"partition sum" 
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W N (E,L) = ^ Js(H N -E) 8{L-Y^nx Pi )DrDp 

(2) 
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where Dr = J\i=i ^ r i; Dp — Ylf^idPi/h 3 ), and e is a 
constant that makes Wn dimensionless. According to 
Boltzmann, the entropy is given by 



S N (E,L) = lnW N {E,L) 



(3) 



Our aim is to find the density profiles that maximize S 
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Following Laliena [Q, the integration over momenta 
in (Q) can be carried out using a Laplace transform. This 
yields 



W N (E, L) - 
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where A is a constant and I = is the inertia ten- 

sor, with elements I a b({r l }) = rnJ2i=i( r i$ab - n^r^b) 
(a, b =1,2, 3). 

In order to evaluate the integral over V , we use the 
following mean field approximation. Letting p(r) denote 
the particles' density inside V (J p(r)dr = N), we set 

*({r,» - *\p] = J f^drdr' (5) 

Iab({ri}) -» Iab[p] = m p(r) {r 2 5 ab - r a r b ) dr (6) 



so that (|J) can be re-cast in the form of the functional- 
integral 



W N (E,L) 
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-P{p]dp(r) 
(7) 



where P[p] is the probability to observe a density profile 
p = p(r). To estimate P[p], we adopt the logic of Lynden- 
Bell 0. We subdivide the spherical volume V into K 
identical cells labeled by the position of their centers. 
The idea is to replace the integral over V with a sum over 
the cells. In order to avoid overlapping we assume that 
each cell may host up to no particles (1 <C no <JV). This 
condition is essentially equivalent to introducing a hard- 
core for each particle, and projects out all the physics 
that is expected to play a role at short distances. P[p] 
is now proportional to the number of ways in which our 
N particles can be distributed inside the K cells with 
maximal capacity no. Denoting by n(r^) the number of 
particles located inside the fc-th cell, a simple combina- 
torial reasoning M leads to 



P[p] oc AH Yl 



cells fc 
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(8) 



Defining the relative cell occupancy c(r) = n(r)/no, so 
that p(r) — Knoc(r) /V, and applying Stirling's formula 
to approximate the factorials (assuming nir^) 3> 1 and 
no — n(rk) ^> 1), one obtains 

P[c] OC Nl e~ ,2 v^ /[c(r) log c(r)+(l-c(r))log(l-c(r))]dr _ 

= Nl e ~^> ^ c( - x ) losc( - x ^ +( - 1 ~ c( - x ^ l ° s( - 1 ~ c( ' a: ^ da: (9) 

where we introduced the dimensionless variable x = r / R 
and defined = n ^ R i = J c(x)dx. 

Plugging this expression into um one easily arrives at 



W N (E,L] 



= S n (E ,L) 



where s% and S2 are given by 
3. 1 
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Sl [c] = ^og[E--L T l- 1 L-^[cl 



(11) 



S2 [ c ] = _ 9 J[c(x)logc(x) + (1 - c(x))log(l - c{x))]dx 

Notice that I = I[c]. We have neglected terms appearing 
in the exponent which do not scale with N. 

For large A~ one can resort to the saddle-point method 
to evaluate (|l0[). Variation of the entropy Sn with re- 
spect to the relative cell occupancy c, with the constraint 
on O enforced by a Lagrange multiplier p, leads to the 
stationarity condition 

log - ~U{x) + xxf-p (12) 

1 — c(xj fc) 2 

or, equivalently, 

C(x) = (1 + e §U( X )-iP(u,x x ) 2 +nyl (13) 

where u> is the angular velocity (related to the total an- 
gular momentum by the relation L = Iw), and the short- 
hands (3 and U (x) are respectively defined as 



= 
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The essence of our mean-field approach becomes clear if 



we notice that $[c 
Eqs ([§) or (£f) 



oc J c(x)U{x)dx. 
are our central result. They hold 
for any long-range potential $ for which our mean-field 
approximation can be justified. Functions c(x) solving 
them and corresponding to entropy maxima represent our 
desired equilibrium distribution of particles. Of course, 
for fixed energy and angular momentum there may exist 
different solutions, each having its entropy. In such cases, 
the higher the entropy the more probable the solution. 

Technically, one can only hope to solve e.g. jl^ ) by 
numerical integration. However, the implicit dependence 
of U (x) on c(x) via a 3-dim. integral makes this a 
formidable task. To simplify things, we pass to spherical 
coordinates, x = (x,9,4>), and expand our potential and 
relative occupancy in series of real spherical harmonics: 
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with x V x 1 = min{a;, x'} and x Ax' — max{i, x'}. bi m {x) 
is a radial function whose precise form we shall soon de- 
rive. Using the series (|l6|), together with the complete- 
ness relation for our basis set {Y; m }, one can easily show 
that U{x) = J2i. (x)Y lm (d,<p), with 



oc / e w ( s iM+M C ]) rfc ( x ) ( 10 ) Ulm ( x ) = 



21 + 1 J (xAx')W lm{ ,[ ' 1 ' 
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Multiplying both sides of (JL3J) by Yi m and integrating over angular variables one finds 



J 



h m (x) = I Y t 



1 _|_ e # ES EL=-i W,0)-^ 2 z 2 sin 2 e+ M 



(18) 
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(a) E = -0.72, L = 0.4 



(b) E = -0.06, L = 0.4 



FIG. 1: Typical behaviour of the radial part bi m as a function 
of I for m — 2. This plot was obtained for E — —0.18 and 
L = 0.44. 
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FIG. 2: Global phase diagram of the self-gravitating JV-body 
system with angular momentum. The dashed lines delimit 
the region where TL[Sn] was computed. 



where I = 0,1, 



and 



= -1,-1 + 1,. 



The 



integral-algebraic system (18) can be solved numerically 
as follows. Fixing E and L and starting from a rea- 
sonable initial guess for 6/ m (a;), one can compute ui m {x) 
from (J!?]) (1-dim. integral). Using this, ( |l8|) can be calcu- 
lated (2-dim. integral) to obtain a better form for bi m (x). 
This scheme can be iterated until convergence. Clearly, 
numerical calculations must be performed with a finite 
number of harmonics. A first reduction is obtained by ex- 
cluding odd harmonics from the calculation. Exclusion of 
/ = 1 fixes the center of mass in the origin, while absence 
of higher-order odd harmonics prevents the formation of 
asymmetric structures (e.g. two clusters of different sizes 
lying at different distances from the origin). Besides this, 




(c) E = -0.9, L = 0.5 




(d) E: 



-0.42, L = 0.5 



FIG. 3: Examples of equilibrium distributions c(x) occurring 
inside our spherical box. Shown are the contour plot and, 
above it, the density profile. The top row corresponds to low 
angular momentum: one sees single clusters (lower E, (a)) 
and disk-like structures (higher E, (b)). The bottom row 
displays some results for high angular momentum: double 
clusters (lower E, (c)) and rings (higher E, (d)). L lies along 
the vertical axis. 



from the typical behaviour of bi m (x) obtained from the 
above procedure, shown in Fig. |l|, one clearly sees that 
bim dies out fast as I increases. Therefore, all calculations 
were performed with even harmonics up to and including 
Z = 16. We set the particles' masses to one and mea- 
sured energy and angular momentum in units of GN 2 /R 
and (RGN 3 ) 1 / 2 , respectively. Solutions of @ at fixed 
E and L were obtained, choosing 9 = 0.02 always. For 
the sake of simplicity we took the angular momentum 
parallel to the 3-axis. 

The resulting phase diagram is reported in Fig. p[ In 
order to discern pure phases from phase-coexistence re- 
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FIG. 4: Entropy as a function of the order parameter In —I22 
at L — 0.5 and different values of E. The values of E and L 
correspond to the four markers ( x ) shown in Fig. ^| 



gions (mixed phases), the Hessian of Sjy, namely 



H[S N ] =det (g E d S a? c 



(19) 



must be analyzed [Q, 01. Indeed, in pure phases one has 
HISn] > 0, while in mixed phases H.{Sn] < 0. In the lat- 
ter, the specific heat is negative and statistical ensembles 
(microcanonical and canonical) are inequivalent. For low 
angular momenta, the system is more likely to be found 
in a single dense cluster (SC) at low energies, while for 
higher energies the most probable state is a gaseous cloud 
(G). The two phases, both pure, are separated by a mixed 
phase with negative specific heat where different equally- 
probable equilibrium configurations coexist. One thus 
recovers the usual collapse scenario Q that is found in 
theories without angular momentum. For higher angular 
momenta, instead, the most probable equilibrium con- 
figuration is a double cluster (DC, pure phase), although 
the gas remains the most probable at sufficiently high en- 
ergies. A sample of equilibrium density profiles is shown 
in Fig. ||. A central question is clearly that of stability 
of such structures. A detailed analysis will be given else- 
where, however rings turn out to be unstable, at odds 
with single and double clusters. 

The appearance of the previously-unobserved double- 



cluster solutions at high angular momenta is particularly 
remarkable (and stresses again the importance of angular 
momentum in self-gravitating systems) , because in such a 
state rotational symmetry is spontaneously broken. This 
is shown explicitly in Fig. ^, where the entropy is plotted 
as a function of the order parameter In — I22 (measured 
in units of NR 2 ), I a b being the components of the inertia 
tensor I (a, b = 1,2,3). (We remind the reader that we 
chose the angular momentum to be parallel to the 3-axis.) 
If L = (or in absence of angular momentum) the system 
is isotropic (In = I22 = I33) and rotational symmetry 
cannot be broken. When L ^ 0, anisotropics may oc- 
cur (J33 ^ I\ 1,222) and one can have either rotationally- 
homogeneous (In = I22) or rotationally-heterogeneous 
(In ^ I22) solutions. The latter correspond to double 
clusters. Indeed, the entropy profile develops two peaks 
at non-zero values of In — I221 corresponding to binary- 
star systems, with the two stars either aligned on the 
1-axis or on the 2-axis. 

Summarizing, we calculated the static equilibrium den- 
sity profiles of a self-gravitating and rotating gas us- 
ing a microcanonical mean-field approach, showing that 
the formation of double-cluster structures (previously ob- 
served only through dynamical approaches) can be ob- 
tained from the spontaneous breakdown of a fundamen- 
tal symmetry of the Hamiltonian (Q). We have presented 
for the first time the global phase diagram as a func- 
tion of the conserved quantities E and L. The inclusion 
of angular momentum in this analysis is absolutely cru- 
cial for these results, the formation of a double-cluster 
structure being possible only through the spontaneous 
breaking of a fundamental symmetry of the Hamilto- 
nian (i.e. the rotational symmetry). However, it goes 
by itself that the formation and stability of stars, dou- 
ble stars, etc. involves other forces (e.g. nuclear and 
sub-nuclear) than Newtonian gravity, and hence could 
require other ingredients than just energy and angular 
momentum juj . The results presented here provide nev- 
ertheless the most detailed static analysis of the problem 
to date, and bridge the existing gap between static and 
dynamical theories. A more exact picture of the situation 
could be obtained by introducing correlations, which are 
ignored in the present mean-field approach. 
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